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Using a symmetry-based atomic scale theory of lattice distortions, we demonstrate that elastic 
textures, such as twin and antiphase boundaries, can generate intricate electronic inhomogeneities 
in materials with strong electron-lattice coupling, as observed in perovskite manganites and other 
functional materials. 
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Recent advances in imaging techniques have re- 
vealed the presence of rich elastic textures in func- 
tional materials such as colossal magnetoresistive (CMR) 
manganites^A^ ferroelectrics,^ ferroelastics,^ and shape 
memory alloys.^ In particular, recent experiments on cer- 
tain perovskite manganites compounds^ have shown the 
correlation between electronic transport properties and 
the presence of meandering antiphase boundaries (APBs) 
within insulating charge ordered domains, interpreted as 
the existence of metallic regions forming around APBs. It 
is also reported that strains near grain boundaries in thin 
film can modify electronic properties in manganites A The 
interplay between elastic texture and electronic hetero- 
geneity is thus central to understanding multiphase co- 
existence and the resultant electronic properties in CMR 
and other functional materials. 

In this Report, we illustrate the importance of elas- 
tic inhomogeneities in modifying electronic properties in 
materials with strong electron-lattice coupling. In par- 
ticular, we study the electronic properties of APBs and 
twin boundaries (TBs) on a two-dimensional (2D) lattice. 
We first show how a recently-developed symmetry-based 
atomic-scale theory of lattice distortions can be used to 
find atomic configurations of twin and antiphase bound- 
aries. Within our framework, we illustrate the differ- 
ences and similarities between TBs and APBs from the 
point of view of localization of long [short] wavelength 
modes inside APBs [TBs], evolution upon energy relaix- 
tion and roughness^ [smoothness] of APBs [TBs]. We 
subsequently perform a tight binding calculation with a 
Su-Schrieffer-Heeger (SSH) type model of electron-lattice 
coupling to predict the distribution of electronic density 
of states, which can be related to the results of scanning 
tunneling microscope (STM) measurements. Our work 
thus forms the basis for predicting electronic properties 
from predesigned materials microstructures. 

Twin boundaries separate domains related by the rota- 
tion of crystalline axes, whereas APBs represent bound- 
aries at which the sequence of alternating distortions, 
such as alternating rotational directions of oxygen oc- 
tahedra in perovskite oxides, change their phase (i.e., 
broken translational symmetry). Although our method 
can be applied to 2D or 3D lattices with mono- or mul- 
tiatomic bases, we illustrate our ideas with a square 
lattice in 2D space with a monatomic basis, for which 
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FIG. 1: Distortion modes for a square lattice in 2D with a 
monatomic basis. ^ For example, Sx(%) — [(i^(i) — (i'^(i + (10)) + 
d^(?+ (11)) - (Tii^ (01))]/2, where d^(?) and (1^(1) represent 
the displacement of the atom at site % along x and y directions 
respectively. 



the appropriate atomic scale distortion variables are the 
modes shown in Fig. Q (Ref. 8). These modes have im- 
portant advantages over displacement variables - they 
reflect the symmetries of the lattice and can serve as 
order parameters (OP) in structural phase transitions. 
Therefore, energy expressions with desired ground states 
can be written in a simpler way in these variables than in 
usual displacement variables. Moreover, since the lattice 
distortions are decomposed into the modes at ^ = (0,0) 
(long- wavelength or intercell modes) and k = (tt, tt) 
(short- wavelength or intracell modes), the approach us- 
ing these modes shows the differences between long and 
short wavelength lattice distortions in a natural way. 

We consider APBs, such as the one shown in Fig. [21 
where open circles represent the distorted atomic posi- 
tions, for Sx or Sy modes. The simplest energy expression 
yielding a ground state with either pure or Sy mode 
lattice distortion is 

i 

+ ^S,(jfSy(if+ Y: (1) 

where the coefficients Gi, G2 (with Gi < G2) can be 
obtained from interatomic forces and Cn are the associ- 
ated elastic moduli. 

Since the ground state has k = (tt, tt) com- 
ponent lattice distortions, we define variables 
with tilde by multiplying (— i)*^+*y with all 
intercell/intracell/displacement variables [e.g., 

5:c(?)=5x(?)(-l)*-+'^ dx{i)=dx{i){-iy-^H In 

wave vector space, this corresponds to interchanging the 



2 




FIG. 2: (Color) An example of the atomic displacement pat- 
tern on either side of an APB. Solid and open circles repre- 
sent the atomic positions for undistorted square lattice and 
distorted lattice, respectively. Red and blue represent the 
positive and negative signs of the Sx mode defined for each 
plaquette. 

(0, 0) and (tt, tt) points. Therefore, Sx and Sy are the 
modes near q = (0,0), and ei, 62, and 63 are the modes 
near q = (7r,7r), where q represents the wave vector for 
the variables with tilde. The constraint equations^ in 
the continuum limit are ei,3(r) = (VySx ± \/xSy)/{2\^) 
and 62 (r) = {VxSx-\-^ySy)/{2^/2). These relations show 
that spatial variations of intracell modes {sx^ Sy) always 
generate intercell modes (ei, 62, 63), which becomes 
microscopic origin of the phenomenological gradient 
energy term in continuum elasticity theory. Therefore, 
intercell modes are present inside an APB at which 
an intracell mode with tilde changes its sign. Similar 
arguments apply to a TE^ - the spatial variations of 
intercell modes always generate intracell modes inside 
TB through their constraints - and lead to a similar 
solitary profile for TB and APB AS However, there exists 
a fundamental difference between intercell and intracell 
modes: In the k ^ limit, ei, 62, and 63 are given as a 
first derivative of dx and dy^ whereas in the q ^ limit, 
Sx and Sy are 2dx and 2dy. Therefore, ei, 62, and 63 are 
related by constraint equations, whereas Sx and Sy are 
not constrained by each other. 

We consider the APB be tween the two ground states 
with Sx{i) = ^^JD IG\ (Fig. [2|). On expanding 
Eq. (QJ around the ground state and retaining leading 
order terms, we write the energy £^apb = ^apb,op + 
^APB,NOP, where 

^;APB,N0P = §^e~n(^)' + E Y^^^^^' 

i,n=l,2,3 % 

i?APB,op = E + (3) 

i 

and D' = D{G2/Gi - 1). The energy ^apb has the 
identical form to that for a TB with 63 as OP^ - a 
double-well potential for OP and harmonic potentials 
for non-OP. For both TB and APB problems, the non- 
OP energy terms mediate an anisotropic interaction be- 
tween OPs through the constraints. The minimiza- 
tion of £^APB,NOP using Lagrange multipliers leads to 
^APB,NOP = T^q^^x{-q)U{q)Sx{q)^ where U{q) has a 
q^ leading order term with an anisotropic coefficient 
U2{0q) = [{Ci + C2 + Cs) + (C2 - Ci - Cs) cos2^,]/16. 
When transformed into real space, the q^ leading order 
term gives rise to a short range interaction between 



OP Sx, where R is the distance between two sites, unlike 
the long-range interaction between OP 63 for the 
TB case.^^'^^ Such different range of interaction is con- 
sistent with well-known understanding^ of smoothness of 
TB and roughness of APB, which our simulation repro- 
duces below. 

The physical origin of the long-range interaction be- 
tween intercell modes and short-range interaction be- 
tween intracell modes lies in the difference in the symme- 
try operations relating the two domains separated by a 
TB or APB. For a TB, the sign change in OP corresponds 
to a change in orientation, which has no intrinsic length 
scale and thus gives rise to a long range interaction. For 
an APB, the sign change in OP signifies translation of the 
configuration by one atomic spacing, implying the pres- 
ence of an intrinsic atomic length scale that is responsible 
for the fast decay of the interaction between OPs. These 
considerations apply to any 2D/3D lattice with mono- or 
multiatomic bases. We can also expect it to be valid for 
the rough and fluctuating stripes suggested in high-Tc 
cupurates, which are examples of magnetic APBs. Ar- 
guments similar to those in Ref. 12, related to the range 
of the interactions, explain the characteristic smoothness 
and roughness of TBs and APBs. 

The full expression of the kernel U (q) in £^apb nop 
given by 

m = ^^Ty{q)-^^Tx{q)^ (4) 

_ {Cr^C2-C^?Tx{q)^Ty{q)^ 

AD' + 2(Ci + C3)T,(g-)2 + 2C2Ty{q)^ ' 

where Tx{q) = tan(gaj/2) and Ty{q) = tdin{qy/2). We 
choose similar parameter values for both APB and TB 
cases and use the Euler method^^ to relax the lattice 
starting from random initial conditions to see difference 
in evolution between TB and APB. The results, 63 for TB 
and Sx for APB, on a 64 x 64 lattice with periodic bound- 
ary conditions are shown in Figs. Ol and 21 respectively, 
where the main panels show real space distributions and 
insets the k or g-space distributions. Figures El^a) and 
Efa) correspond to early stages of the relaxation, which 
show the characteristics of the OP distribution at high 
temperature(T), whereas Figs. |3fb) and^b) reflect the 
OP distributions for late stages at low T. Even at the 
early stage, the presence of the long range correlation be- 
tween intercell OP along 45^ and 135^ can be identified 
in the main panel in Fig.[3fa), reminiscent of tweed struc- 
tures in martensitic materials. The /c-space distribution 
plotted in the inset of Fig.|3Ja) shows strong preference of 
es{k) along 45^ and 135^ orientations. Such long-range 
correlation and anisotropy are absent for the phase of 
the intracell mode distortion, as shown in Fig.^a). The 
late stage of the relaxation depicted in Figs. EJb) and 
^b) show characteristic smooth TBs^ and rough APBA 
The TB is metastable and cannot be removed by fur- 
ther relaxation, unless large noise is applied. In contrast, 
the APB in Fig. ^b), which is a ring due to the peri- 
odic boundary condition, shrinks and disappears upon 



FIG. 3: (Color) 63(1) (main panels) and |e3(/c)| (insets) ob- 
tained from simulation for TB: (a) early and (b) late stages 
of the lattice relaxation. Parameter values (see Ref. for the 
definitions) are Ai = = B = 4, A3 = 5, and F3 50. The 
center of the inset corresponds to ^ = and the four corners 
k = (±7r,±7r). (Dark red: 0.32, dark blue: -0.32, white color 
in the inset: larger than 0.2 times the maximum of |e3(A:)|). 



FIG. 4: (Color) Sx{i) (main panels) and |scc(^| (insets) dis- 
tributions obtained from simulations for APB: (a) early and 
(b) late stage of the lattice relaxation. Parameter values are 
Ci = C2 = C3 = = 4, D = 5, and Gi = 50. 



further relaxation. Although the solitary-wave profile of 
the smooth APB along a certain direction is a metastable 
state, ^ the absence of a long range interaction between 
the phase of intracell mode distortions prevents the re- 
laxation of random initial configuration from reaching 
such metastable states. This indicates that lattice de- 
fects or boundary conditions are necessary to reach the 
metastable solitary- wave APB configuration in materials, 
and influence the geometry of APBs. 

To study functional electronic aspects associated with 
these TB and APB microstructures, we consider the 
modulation of electronic properties based on the follow- 
ing SSH electron-lattice coupling Hamiltonian, 

i^ssH= J2 -*o[l-a(rff+,-d?)](4c.+a+4+a^?)' (5) 

i,a=x,y 

where cl is the creation operator for an electron at i. We 

use d^{i) and d'^{i) obtained from our atomistic model as 
inputs to the SSH Hamiltonian. For to = 1, a = 1, and 
the TB and APB results shown in Figs. [3Jb) and^b), 
we find all energy levels and eigenstates numerically, and 



FIG. 5: (Color) Electronic properties calculated for the lat- 
tice distortion in Fig. I^b). (a) Inset: local DOS within the 
domain (blue curve) and TB (red curve); main panel: spatial 
distribution of local DOS at Ep. (b) Corresponding charge 
density distribution (dark blue: 0.36, green: 0.39). 




FIG. 6: (Color) Electronic properties calculated for the lat- 
tice distortion in Fig. 2Jb). (a) Inset: local DOS within the 
domain (blue curve) and APB (red and green curves); main 
panel: spatial distribution of local DOS at Ep = (dark 
blue: 0, dark red: 2). (b) Spatial charge density distribution 
for Ef = 0.1 (green: 0.52, dark red: 0.72). 



calculate the local density of state (DOS) at each site 
and the distributions of local DOS and charge density 
for chosen Fermi energies {Ep)^ as shown in Figs. El 
andlH 

The blue curve in the inset of Fig. [3ja) represents the 
local DOS within the domain of Fig. IHtb), close to the 
bulk DOS for the homogeneous phase. The energy dif- 
ference between the two peaks is proportional to a|e3|. 
Within the TB, \e^\ is small and therefore the local DOS 
spectral weight moves toward £^ = 0, as shown with a 
red curve. This local shift of the DOS weight can be 
measured using direct local probes, such as STM.^'^^ For 
Ep shown in this inset, the local DOS is smaller inside 
the TB than inside the domain, and the main panel of 
Fig. El^a) shows the real space distribution of local DOS 
at this Ep- The oscillation of local DOS within the do- 
main is related to Friedel oscillations. The feedback from 
electron to lattice neglected here may generate a simi- 
lar Friedel oscillation of the lattice. The charge density, 
which is the area in local DOS below Ep^ is a constant 
(0.5) if Ep = 0. If Ep < 0, more electrons are depleted 
from the TB than the domain, as shown in Fig. ISfb). 
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The charge density also has a Friedel oscillation, though 
not clearly visible in Fig.jSfb). 

The inset of Fig. Efa) shows the local DOS for the 
APB case. Within the domain, each site has a neighbor 
along X axis that is closer in one direction and farther in 
the other (Fig. (2)). Such a distortion pattern results in 
a V-shaped local DOS with zero DOS at £^ = 0, plotted 
as a blue curve in the inset of Fig. Efa) (the small finite 
value of local DOS at £^ = is due to the finite energy 
bin size). Within the APB, the average of the two bond 
lengths along x-axis has a (tt, tt) component. For exam- 
ple, along sites at the center of the 90^ APB shown in 
Fig. El the neighboring sites in the x-direction are alter- 
nately closer or farther apart. The red line in the inset 
of Fig. ISfa) represents the local DOS for the sites with 
two farther neighbors within the center of APB. This has 
a peak near E = {). The green line represents the local 
DOS at sites with two closer neighbors, right next along 
y direction to the sites for the red curve. The difference 
between the red and green curves emphasizes the differ- 
ence that occurs in density of states due to atomic-scale 
changes in microstructure as along an APB. Such atom- 
istically sharp changes at structural interfaces have been 
seen in CMR manganites with STM (Renner et al}). 
The local DOS at = is zero within the domain 
and finite only around the APB with a (tt, tt) component 
modulation, as depicted in the main panel of Fig. [Sfa). 
This means that the electronic states created by the APB 
dominate the low energy properties of the whole system, 
e.g., conductivity or specific heat. The charge density 
for Ep slightly higher than zero is plotted in Fig. IHl^b). 
Checkerboard-type modulations of electronic properties 
around APBs, refiecting the electronic states coupled to 
the (tt, tt) type lattice distortions, demonstrate that STM 
studies of TB and APB can be used to reveal underly- 
ing electronic properties of materials (usually hidden in a 
homogeneous phase), similar to the way single impurity 
atoms have been used to reveal superconducting proper- 
ties in high Tc cupurates. Results for different values 
of Ep are displayed in Fig. [3 Local DOS modulations, 
particularly the checkerboard type modulations, tend to 
reach farther into the domain than the elastic texture 
itself. We also note that the checkerboard type modu- 
lations with a wave vector (tt, tt) are present at any Ep^ 
in contrast to Friedel oscillations for TBs which change 
wavevectors with Ep- 

Although the orbital states and type of electron-lattice 



coupling in CMR manganites are quite different from the 
simple model presented here, our results suggest that 
simple APBs from double-well type potentials, such as 
Eq.(3), may give rise to electronic heterogeneities, but 
not undistorted metallic regions with uniform charge den- 
sities near APBs. An energy landscape with a local min- 
imum at undistorted state, such as the one considered 
in Ref. 14, would be more likely to nucleate metallic do- 
mains at APBs, and create percolating conducting paths 
in CMR manganites. 



In summary, we have shown that in functional mate- 




FIG. 7: (Color) Local DOS calculated for the lattice distor- 
tion in Fig. 2Ib) at various Bp values (a) -1.5, (b) -1.1, (c) 
-0.7, and (d) -0.3. (dark red: 0.35, dark blue: 0) 



rials with strong electron lattice coupling, the electronic 
properties are modified near elastic textures such as TBs 
and APBs, which can be directly measured by STM. The 
results also show that the heterogeneities of electron lo- 
cal DOS are not just confined within TBs and APBs, but 
can propagate into domains in the form of Friedel oscil- 
lations for TBs and with the wave vector of short wave 
length lattice distortions for APBs. 
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